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Abstract 

The pion number fluctuations are considered in the system of pions and large mass fireballs decaying 
finally into pions. A formulation which gives an extension of the model of independent sources is 
suggested. The grand canonical and micro-canonical ensemble formulations of the pion-fireball gas 
are considered as particular examples. 
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1. During several decades the models of hadron productions in high energy collisions from 
decays of the fireballs I], droplets 2], strings [3], quark-gluon bags [4] have been considered. 
We will use the name fireball for all these objects with high mass m. A typical estimate is 
m > 2 GeV. Fireballs may crucially influence the thermodynamical properties of the system. 
For example, their presence leads to the limiting temperature due to the exponential mass 
spectrum of fireballs [1] and strings [3], or to the phase transitions in the gas of quark-gluon bags 
[4j. It is assumed that the fireballs decay into stable hadrons, mainly to pions. A comparison 
with the data on hadron multiplicities at high energy collisions does not indicate a presence 
of massive fireballs (see however recent discussion in Ref. [5[]). Thus their contribution to 
the hadron average multiplicities is probably rather small, if it exists. On the other hand, 
even for small contributions of fireballs to average pion multiplicity, there is a chance to find 
the signatures of large mass fireballs by measuring the multiplicity fluctuations. This will 
be discussed in the present paper. The study of event-by-event fluctuations in high energy 
nucleus-nucleus collisions may also open new possibilities to investigate the phase transition 
between hadronic and partonic matter as well as the QCD critical point (see, e.g., the reviews 

The scaled variance 

(N 2 ) - (N) 2 



(N) ' 

is a useful measure to describe the pion multiplicity fluctuations. In Eq. (pQ), N is the total pion 
number after fireballs decay. The event-by-event averaging is defined by the pion probability 
distribution W(N), 

oo 

(...) = j\..W(JV). (2) 

N=0 

If there is only one type of fireballs, and they are the only sources of pions, the scaled 
variance ui is equal to (see, e.g., 



u f + n f u* f , (3) 



where n/ and ujj is, respectively, the average number of pions and the scaled variance of 
pion multiplicity fluctuations from one-fireball decay. In Eq. ([3]), uj*j is the scaled variance 
for the fluctuations of the number of fireballs. The Eq. ([3]) corresponds to the contributions 
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of independent sources (fireballs). In this paper we present the extension of the model of 
independent sources in two directions. First, we consider both fireballs and primordial pions 
on an equal footing. Thus, only a part of the final pion multiplicity (probably a small part) is 
due to the fireball decays. Another part is just the number of the primordial pions. Second, 
only the decays of fireballs (sources) are assumed to be independent, while their production 
may include different type of correlations. These correlations are absent in the standard model 
of independent sources. 

2. We introduce the probability distribution function P(N , N%, . . . , N k ) to describe the 
production of fireballs and primordial pions. The distribution function Tj(n) describes the 
decays of i-th fireball into pions. The probability W(N) to find N pions in the system of 
primordial pions and fireballs decaying finally into pions is then calculated as: 

oo k Ni oo / k Ni \ 

w(n) = Wb^,...,iv fc ) niIE r ^) 5 p-EE^ • (4) 

N ,Ni,...,N k =Q »=0 3=1 n»j=0 \ i=0 j=l / 

In Eq. (JIJ), i = corresponds to pions and i = 1, k to different types of fireballs, N is the 
number of i-th species, and is the number of pions from the decay of fireball of type i with 
the number j = 1, . . . , jVj. The normalization conditions, 



P(N ,N u ...,N k ) = 1, ^r,(n) = l, (5) 

N ,Ni,...,N k =0 n=0 

are assumed, and they lead to J2n W(N) = 1. The following assignment will be used: 

oo 

(...)t = E ■■■P(N ,N 1 ,...,N h ) . (6) 

N ,Ni,...,N k =0 

According to (Jl]) the fireballs decay independently, i.e. the number of pions from the decays 
of two different fireballs do not correlate. One then finds, 



oo 



2j Yi{nij)nij = rii , ^2 ^i( n ij)^p( n pk)nijn pk = rnn p + 5 ip 5 jq utrii , (7) 



ny=0 iiij,n pk =0 

where 



1 71=0 
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with no — 1 and loq — 0. For convenience the pions are considered as fireballs decaying into 
themselves, i.e. To(riij) = 5(riij — 1). Introducing the average numbers, (Ni) T = Ni, and 
correlators, (AA^AA(j} T , where AN = N^ — (Ni) T , one obtains: 

(AO = ^iV^, <(AiV) 2 ) = ^u i N i n i + Y,n i n j (AN i AN j ) T , (9) 

i=0 i=l i,j=0 

where A A" = N — (N). In Eq. ([9]) the terms uJiNiHi describe the pion fluctuations at fixed 
number of fireballs Ni due to the probabilistic character of fireball decays. The next terms 
for ((AN) 2 ) in Eq. (Q present the contribution due to the fluctuations and correlations of N 
numbers. Using Eq. 0, the scaled variance ([T]) is presented as: 

u = ^fc - w • ( 10 ) 

E i= o n i N i 

The same structure of equation for u was previously obtained for the equilibrium hadron- 
resonance gas within the canonical ensemble jl0|] and micro-canonical ensemble Uj. The form 
of Eqs. ([9]) and ( flOl) appear to be rather general and their validity do not require any thermal 
equilibrium, i.e. the averaging (. . .) T may include arbitrary dynamical (non-thermal) effects. 

3. If no fireballs exist, Eq. ( TlOl) is reduced to 

w = = = • (11) 

1\ o 



If there is only one type of fireballs and no primordial pions, i.e. k = 1 and A^o = 0, Eq. (TTOT) 
reads: 

_ ((AiV 1 ) 2 ) r _ , 

U = U>\ + Tlx = = Ui + Hi Wj , (12) 

and corresponds to the model of independent sources (J3j). 

Let us consider now the system when both the primordial pions and fireballs are present. 
For one type of fireballs and no correlations between fireballs and primordial pions, i.e. 
(AA r AA ? " 1 ) :r = 0, the scaled variance ffTUj) reads: 

u _ 1+ K - I) N Q +^tNinj_ + gigi (ui - 1) 

N + JV x ni 

The Eq. ( fl3l) is reduced to ( TTTI) and ( |T2l) at ni = and A^ = 0, respectively. Let us fix the 
ratio R of 7r-mesons from the fireball decays to the total number of 7r-mesons, 

R = _ . (14) 

N + n 1 N 1 
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Using Eq. ( JT4l) . the scaled variance (1T31) can be rewritten as, 

u = Wq{1-R) + (ooi + niul)R. (15) 

The Eq. (115p includes the results for the pure pion system ( TTTT) and pure fireball system ( [121) 
as the special cases at R = and -R = 1 respectively. Assuming 

ni(m) = 7 m , (16) 

where 7 _1 has a meaning of the average pion energy from the decay of a fireball, one finds 
that u increases linearly with fireball mass. This is shown in Fig. [H where the parameters 
in Eqs. (I15II16P are fixed as uo^ = uo\ = uj\ = 1 and 7 _1 = 0.5 GeV. If all pions come from 
the fireball decays, i.e. R — 1, the scaled variance u becomes more than 10 times larger than 
Uq — 1 at m = 5 GeV. Even for very small average pion multiplicity from fireball decays, i.e. 
when R 1, the fireball contribution to u always dominates at large m. 




1 ■ ' ■ ' ■ ' ■ 1 

5 10 15 20 
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FIG. 1: The scaled variance (115[) as the function on m is shown at R = 0, 0.05, 0.1, 0.25, 0.5, 
and 1. The lower bound for the fireball mass is taken as m = 2 GeV. The parameters are fixed as 
= u* = loi = 1 and 7 _1 = 0.5 GeV. 
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4. The statistical mechanics of pion-fireball gas in the grand canonical ensemble (GCE) 
gives an example when (ANqANi) t vanishes. This leads to Eq. ( fl3l) . In the Boltzmann 
approximation the distributions of N and Ni in the GCE are the Poissonian ones, hence 

= u* = 1, like they are fixed in Fig. [TJ The average particle numbers N Q and Ni are: 



W„(V,T) = Y>!efK,(!SZ) *%VT», (17) 
N l( V,T) = ^L Kl (Z) a 9l C ^V /2 exp f- ^ , (18) 



2tt 2 VT/ 3 V 27r / V T - 

where V and T are the system volume and temperature, m , go and m, pi are the mass and 
the degeneracy factor of the pion and fireball respectively. The energies E , E\ and specific 
heats C , G\ of the pions and fireballs are respectively, 

dN — - - dN l ( 3 \ — 



E = T 2 ^ = 3TN , Ei = T -^f = m + ^T)^, (19) 



C = ^ - 12iV , Gl ^ 



m\ 2 m 15 

r) +3 t + T 



iVi . (20) 



dT u ' 1 <9T 

The non-relativistic approximation m/T ^> 1 for fireball thermodynamical functions and ultra- 
relativistic approximations itlq = for pions are presented at the last steps in Eqs. ( 11711201) . 
They will be used further to obtain simple analytical estimates. The thermodynamical functions 
(I18ti20p of the fireball component include the exponential suppression factor exp (—m/T) at 
m/T>l and thus go to zero at large mass. To have the finite values of R at m/T ^> 1 one 
needs large fireball effective degeneracy factor g-y. Using Eq. (117118)) . from Eq. (|14p at fixed 
R one finds g\ ~ m~ 5//2 exp(m/T) at m/T 1. This behavior resembles the fireball mass 
spectrum p{m) = Cm a exp(m/Tn) assumed in the Hagedorn model PJ. The integration over 
fireball masses with spectrum p(m) leads to the limiting temperature T = T# in the Hagedorn 
model. Our model consideration in the present paper is rather different. The fireball mass m 
is assumed to be fixed. There is no integration over the fireball mass spectrum and thus there 
is no limiting temperature in the system. 

For thermal massless pions the average energy per particle equals to 3T. The typical temper- 
ature of the hadron gas is about 160 MeV. Thus the average energy 7 _1 =0.5 GeV of pion from 
the fireball decay assumed in Fig. [1] is close to the average thermal energy of the primordial pi- 
ons. However, other relations are also possible. The resonance decays in the hadron-resonance 
thermal gas at different T give the examples when the energy of pions from resonance decays 
can be both larger and smaller than the thermal energy of the primordial pion. 



The GCE formulation of the pion-fireball gas leads to a linear increase of to with m for fixed 
value of R. Thus large values of u are possible even at small R if fireball mass is large enough. 
This is connected with unusual thermodynamic properties of this system. At finite values of 
R the ratio E\/Eq is also finite. It is proportional to R at R <C 1 and thus is close to zero. 
On the other hand, Ci/Cq is large and goes to infinity at large m and fixed R. Thus fireball 
contribution to the total energy can be rather small (i.e. E\/Eq <C 1 when R <C 1), but the 
specific heat of fireball component always dominates, C\/Cq 3> 1, at large m. 

5. In the GCE, the average energy is fixed and thermal energy fluctuations around this 
average value are admitted. In the micro-canonical ensemble (MCE) 1 the total energy is fixed 
for each microscopic state of the pion-fireball gas. This exact energy conservation generates 



specific corre 
presented as 



ations between particle numbers. The MCE correlators (AiVjAiVj) can be 
(ANAN,)^ = Sij £ < ( - E ^-' j>^" , (21) 



p,l w p,l E p,l 



where e pi = a/p 2 + mf and v pi = n P}i (l ± n pi ) with plus sign for bosons and minus sign 
for fermions. The values H-p j £1X6 the average occupation numbers of states, labeled by 3- 
momentum p. In the Boltzmann approximation the above expressions are simplified: v i pi = 
n P) j = exp(— e P) i/T). In the thermodynamic limit of large volume V, the sums Yl p ■ ■ ■ m Eq. ( 12T1) 
can be replaced by the integrals and presented in terms of the GCE thermo dynamical functions 
([T71I20D: 



= fSf ^ 6XP = ~Ni(V, T) , (22) 



X>p>p, = f£ jf °P 2d P \jP 2 + m * ex P ( ~ V " T - ) = Ei(V,T) , (23) 
= ^[p 2d P (P 2 + ™ 2 ) exp(-^^) = T* Cl (V,T). (24) 




p 

Using Eqs. (122H24P . the correlators ( T2TT) can be presented as 



(AiV.AiV,)^ = SyNi - 7^C i ■ (25) 



1 It can be also named the grand micro-canonical ensemble [12[ as no conserved charges are introduced. 
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The second term in the r.h.s. of Eq. ( 1251) comes due to the exact energy conservation in the 
MCE. It is absent in the GCE, (ANiANj) M = SyNi. Note that being different from the 
GCE values the MCE correlators (|25|) can be expressed in terms of the GCE quantities. The 
temperature parameter T is fixed by the requirement that GCE average energy is equal to 
the fixed value of the MCE energy. Under this condition the average values iVj and are 
the same in the MCE and GCE. The average final pion multiplicity (N) and the pion number 
fluctuations expressed in {(AN) 2 ) or u are given by general formulae (jHJ) and (fTOf) . In Eqs. (19|ll0p 
the correlators (ANiANj) T should be substituted by the MCE ones (125]) . 

Let us start with an example of the pure pion gas with no fireballs. The scaled variance u 
then reads: 



-=2 



a = (WL = t _ , ( 26) 

JV PQN„ 



instead of u = 1 for the Boltzmann pion gas in the GCE. The second term in r.h.s. of Eq. 
leads to a suppression of the pion number fluctuations due to the exact energy conservation in 
the MCE. For massless pions, Eq. (1261) gives uj = 1/4. This result was obtained for the first 
time in Ref. [3]. If there are only fireballs with mass m and no primordial pions (i.e. Nq = 0), 
using Eq. ffTUl) . in the MCE one finds: 

3 /T x 2 




= - m ( 1 - T2 ^— + n x u}{ = u\ + jm - y— J =£ ^, . (2, ) 



where the relation n\ = 7m f fT6|) and non-relativistic approximations for Ni, Ei, and C\ ( 1181 - 
2HD are used. The Eqs.([2SD and (EB give the MCE realization of Eqs. (EETJ and (PJ) respectively. 
The scaled variance for the fireball number fluctuations in the MCE is coC = 3(T/m) 2 /2. This 

n 

was obtained for the first time in Ref. |13[]. For m/T — >• 00 the scaled variance a>* goes to zero as 
m -2 in the MCE. Thus niul goes to zero as m -1 , and the pion scaled variance uj approaches uo\ 
at large m, i.e. it is defined by the only fluctuations due to the one-fireball decay. This behavior 
is completely different from the GCE result. In the GCE, uj* = 1 and the term n\UJ* increases 
linearly with m. Thus the contribution due to the fireball number fluctuations dominates in 
GCE at large m. 

For the case when both the primordial pions and one type of fireballs of mass m exist one 
finds, using Eq. ( fl4l : 

6 7 T(1 - R)R{2jT - 1) + f 7 T(R 2 + 6R-6)- (T/m) 
U = 1 + R{Ul ~ l) + R+[12 1 T(l-R) + 3R]-(T/m) ' (28) 
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where the relation ni = 7m (|16|) . non-relativistic approximations for Ni, Ei, and C\ (I18H20I) 
and ultra-relativistic approximation for No, E , and Cq (I17II19H20"]) of the pion component are 
used. 




m (GeV) m (GeV) 

FIG. 2: The scaled variance to (|28p in the MCE for the pion fireball gas as the function of fireball 
mass m at different values of jT. The lower bound for the fireball mass is taken as m = 2 GeV. 
The solid line corresponds to R = 0.05, dashed line to R = 0.25, and dotted line to R = 0.75. The 
parameters are fixed as T = 160 MeV and u\ = 1. 

The dependence of u on particle mass m is presented in Fig. [2j One finds that u approaches 
to finite value at m/T —>■ 00. At small 7T the limiting values of u is reached at moderate 
values of m, as it is seen in the upper panel of Fig. [3 Most rapidly the asymptotic value of u 
is reached at 7T = 0.25. It requires much larger m at large 7T as the lower panel of Fig. [2] 
demonstrates. The Fig. [2] also shows that the smaller is R the larger values of m are needed to 
reach the limiting values of uj. The energy of pion from fireball decay should be larger than pion 
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mass. This leads to the restriction 7 < l/m T . As the typical temperature is T = 160 MeV, the 
largest physically allowed value of 7T is approximately equal to 1. 

The limiting values of uj for m/T — > 00 at fixed non-zero R are equal to: 



w = 1 + R {uj! - 1) + 67T (1 - R)(2^T - 1) . 



(29) 



The dependence of the scaled variance (1291) on the dimensionless parameters 7T at different R 
and on R at different 7T is presented in Fig. [3] left and ngM respectively. The parameter oj\ is 
fixed in Fig. CD as uj\ — 1. 




6 

5 

S 4 

3' 
2 
1' 
0' 



0.25 



1 1 ■ 

k. 


Y T=0.25; 




Y T=0.5 ; 




7 T=0.75. 




•\_ yT=1 


i- . . 




\ 




► — 



0.5 

R 



0.75 



FIG. 3: The dependence of the scaled variance u on 7T for different values of R = 0.01, 0.25, 0.5, 1 
(left) and on R at 7T = 0.25, 0.5, 0.75, 1 (right). The Eq. (Tj9|) is used with u x = 1. 

The minimum of ( 1291) corresponds to 7T = 1/4, and the minimal value is uj = 1/4 + R(uj\ — 
1/4). If ui < 1/4, the minimal value of uj is for R = 1, and it equals to u; = at uj± = 0. If 
uj\ > 1/4, the minimal value of uj is for R —>■ and it equals = 1/4. 

At 7T 0, it follows from Eq. (1291. a; = 1 + #(^1 - 1). This is shown in Fig. 1 /e/l If 
i? — >• one obtains uj = 1 for all values of tui. In Fig. [3] the value of uji is fixed as U\ — 1. Thus 
u; = 1 at 7T — > for all values of R. 

The simultaneous limit m — > 00 and -R — > leads to the uncertainty. Introducing mR/T = 
a = const one finds at m — » 00: 

a(2oT-l)6oT-9QT 

a + 12 7 T v ; 

Changing a in the last formula from to 00 one obtains the values of uj from 1/4, which is the 

MCE value for the massless pion gas in the Boltzmann approximation, to 1 + 6 , ~fT(2T'-f — 1). 
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FIG. 4: Left: The dependence of the scaled variance u> on R at 7T — > for different values of 
u)\ = 0, 0.25, 1. Right: The dependence of the scaled variance lo (|3U|) on 7T at different a. The lower 
solid line lo = 1/4 corresponds to a = and the upper line 1 + 6 r yT(2T r y — 1) to a = 00. Possible 
values of cj in the simultaneous limit m — > 00 and i? — > lie between these two lines in the dashed 
region. 



This is shown in Fig. H] ng/ii. The fireballs do not contribute to the pion multiplicity at R — > 0. 
However, their presence affects uj. 

6. Note that the numbers of positively N + and negatively iV_ charged pions are usually 
measured. Effects for the fluctuations of N + and N_ due to neutral vector mesons decaying 



into 7r + 7r -pairs have been discussed in Ref. 



141 ] . Fireballs decaying into large numbers of 



7r + 7r -pairs make these effects very strong. Assuming that fireballs have zero electric charge 
one finds for Q = N + - and N ch = N + + iV_ in the GCE: 

<(AQ) 2 ) = (1 - R)(N ch ) , ((AN ch ) 2 ) = [1 + (u ch -l)R + n ch R](N ch ) , (31) 

where R is the ratio of the number of charged pions from the decays of fireballs to the total 
number of charged pions, n c h ~ m and u c h ~ 1 are respectively the average number and scaled 
variance of charged pions from one fireball decay. A presence of fireballs leads to the strong 
enhancement of the N c h fluctuations and to the suppression of the Q fluctuations. As seen from 
Eq. ([SID, the ratio 

<(AQ) 2 > 1 - R 



((AN ch )*) ~ 1 + (cu ch - 1)R + n ch R 



(32) 
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is sensitive to both R 
may be applicable in 
intervals are satisfied 



and n c h and it goes to zero if R is fixed and m — > oo. The GCE analysis 
ligh energy pion production if the following conditions for the rapidity 

3, 



Ytot > Yacc > 1 , (33) 

where Y tot is the total rapidity interval allowed in high energy collision and Y acc is the rapidity 
interval for the accepted pions. The first of these inequalities is needed to relax the energy 
conservation effects and the second inequality ensures that the pion rapidity spreading from 
fireball decay is much smaller than the accepted interval. 

7. The pion number fluctuations have been considered in the system of pions and fireballs 
decaying into pions. Our model formulation gives a generalization of the model of independent 
sources. Both statistical and dynamical effects in the production of primordial pions and 
fireballs can be included within the present scheme. The statistical equilibrium within the grand 
canonical and micro-canonical ensemble formulations for the pion-fireball gas are considered in 
details as particular examples. Important result of these studies is a strong suppression effect 
for the pion number fluctuations in the MCE due to the energy conservation. In the GCE the 
scaled variance of the fireball number fluctuations is uj^ = 1, and this gives a linear increase of 
the scaled variance u of pion number fluctuations with fireball mass m due to a large number of 
pions from one-fireball decay, n\ ~ m. In the MCE the uj* is strongly suppressed at high masses 
m. This leads in the limit m —>■ oo to finite values of lu depending of the model parameters R 
and jT. 

Even for small contributions of fireballs to the average pion multiplicity there is a chance to 
find the signals of large mass fireballs by measuring the event-by-event multiplicity fluctuations. 
If final pions are detected in a small part of the phase space the exact energy conservation 
becomes not so much important. This leads to a strong enhancement of the pion multiplicity 
fluctuations and linear increase of the scaled variance u with fireball mass m. Both conditions 
f lHHj) can be simultaneously satisfied at high collision energies. They are only approximately 
fulfilled in nucleus-nucleus collisions at RHIC, and can be much stringent in the future LHC 
experiments. 
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